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Abstract 

We call Dyson process any process on ensembles of matrices in which the entries 
undergo diffusion. We are interested in the distribution of the eigenvalues (or singular 
values) of such matrices. In the original Dyson process it was the ensemble of n x n 
Hermitian matrices, and the eigenvalues describe n curves. Given sets Xi, . . . ,Xm 
the probability that for each k no curve passes through at time is given by the 
Fredholm determinant of a certain matrix kernel, the extended Hermite kernel. For 
this reason we call this Dyson process the Hermite process. Similarly, when the entries 
of a complex matrix undergo diffusion we call the evolution of its singular values the 
Laguerre process, for which there is a corresponding extended Laguerre kernel. Scaling 
the Hermite process at the edge leads to the Airy process (which was introduced by 
Prahofer and Spohn as the limiting stationary process for a polynuclear growth model) 
and in the bulk to the sine process; scaling the Laguerre process at the edge leads to 
the Bessel process. 

In earlier work the authors found a system of ordinary differential equations with 
independent variable ^ whose solution determined the probabilities 

Pr {A{ti) < a + e, • • • , A{t^) < U + , 

where r — > A[t) denotes the top curve of the Airy process. Our first result is a 
generalization and strengthening of this. We assume that each is a finite union of 
intervals and find a system of partial differential equations, with the end-points of the 
intervals of the X^ as independent variables, whose solution determines the probability 
that for each k no curve passes through X^ at time r^,. Then we find the analogous 
systems for the Hermite process (which is more complicated) and also for the sine 
process. Finally we find an analogous system of PDEs for the Bessel process, which is 
the most difficult. 
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I. Introduction 



Wc call Dyson process any process on ensembles of matrices in which the entries un- 
dergo diffusion. In the original Dyson process [3] it was the ensemble of n x n Hermitian 
matrices if, where the independent coefficients of each matrix H independently execute 
Brownian motion subject to a harmonic restoring force. In one dimension this is the familiar 
Ornstcin-Uhlenbeck (velocity) process. The solution to the forward (Fokker- Planck) equa- 
tion generalizes to the matrix case with the result that the probability density of H at time 
T — T2 corresponding to the initial condition H = H' aX t = ti is a. normalization constant 
depending upon n and q times 

^ TT{H-qH'f\ 



exp 



where q = e'^^~'^^. As Dyson observed, the equilibrium measure as r2 — > oo is the GUE 
measure of random matrix theory. We refer to this particular Dyson process as the Hermite 
process for reasons that will become clear below. 

With initial conditions at time Ti distributed according to the GUE measure, the prob- 
ability that at times r^. (A; = 2, . . . ,m),^ H{Tk) is in an infinitesimal neighborhood of is 
a normalization constant times 

exp (-IV H^) n cxp LTLiELZllzlElzll] dH,... dH^ , (1.1) 

j=2 \ J 

where qj — e^^~'^^+^ . Alternatively, (1.1) can be interpreted as the equilibrium measure for a 
chain of m coupled n x n Hermitian matrices Hk- 

In random matrix theory, and more generally Dyson processes, one is interested in the 
distribution of the eigenvalues (or singular values) of H. It is a classical result of Gaudin 
[6] that the distribution functions for the eigenvalues in GUE are expressible in terms of 
the Predholm determinant of an integral kernel called the Hermite kernel. In the process 
interpretation, the evolution of the eigenvalues can be thought of as consisting of n curves 
parametrized by time. Given Ti < • • • < and subsets of R, the quantity of interest 
is the probability that for all k no curve passes through at time tj,. It follows from 
the work of Eynard and Mehta [4] that this probability is also expressible as the Predholm 
determinant of an extended Hermite kernel, an m x m matrix kernel related to the kernel 
associated with the random matrix ensemble corresponding to the equilibrium distribution.^ 

Here is how it is derived. One first diagonalizes each Hk and then employs the Harish- 
Chandra/Itzykson-Zuber integral (see, e.g. [18]) to integrate out the unitary parts. The result 
is that the induced measure on eigenvalues has a density P{Xii, ■ ■ ■ , Ai„; . . . ; A^i, . . . , Xmn) 
given up to a normalization constant by 



n e n det (e^'^'^'^"^'^) A{X,) A(A„), (1.2) 

k=l k=l 



^We are assuming ti < ■ ■ ■ < Tm- 

^This was described in a lecture by Kurt Johansson [8], who recently communicated to us a sketch of his 
derivation [9] . Matrix kernels, of a different kind, also appear in [1] . 
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where qo — Qm — and A denotes Vandermonde determinant.^ 

In [4] it was shown that for a chain of coupled matrices with probabihty density of 
this type the correlation functions could be expressed as block determinants whose entries 
are matrix kernels evaluated at the various points, generalizing Dyson's expression for the 
correlation functions for a single matrix. As with the case of random matrices, one could 
then get a Fredholm determinant representation for the probability that for each k no curve 
passes through at time r^. In the case at hand the matrix kernel 

L{x,y) = {L,,{x,y))^^^, 

is the extended Hermite kernel and has entries 

n—l 

J2 e*^ ^^'-^J) ipk{x) (Pk{y) if i > j, 

k=0 

(1.3) 

oo 

- E e^^^^-^i) ipk{x) cpkiy) if i < j. 

k=n 

Here ipk are the harmonic oscillator functions e~^^/^pfc(x) where the pk are the normalized 
Hermite polynomials. If K is the operator with matrix kernel (Kij), where 

Kij (x,y)^ Lij (x, y) Xx^ (y) , 

then the probability that for each k no curve passes through Xk at time is equal to 
det (/ — K). In the special case Xk — {^k, oo) this is the probability that the largest 
eigenvalue at time is at most ^k- 

It is natural to consider also the evolution of the singular values of complex matrices. 
This is the Dyson process on the space of p x n complex matrices. (We always take p > n.) 
The analogue of (1.1) here is [5] 

exp (-Tr AlA,) f[ exp {{A, - - q,A,_,)) \ 

j=2 \ ^ J 

After integration over the unitary parts this becomes a normalization constant times 

(2 \ n 
.X ^'^^ S'" n rtet (/. (t^^Tv:^)) 

k=i k=i \ \^~ik+i )) 

n n 

X A(Ai) A(A^) n Aif n c^An ■ ■ ■ ciA^n, (1-5) 

j=l i=l 

where Jq, is the modified Bessel function and a — p — n. (The \ki are the squares of the 
singular values.) This is of the same general form as for the Hermite process, and here also 
there is a corresponding matrix kernel, the extended Laguerre kernel. It is given by the same 

^This expression shows the connection with the theory of determinantal processes, in which probabihty 
densities are defined by products of determinants [5, 7, 12]. 



Lij{x,y) = < 
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formulas (1.3) as before, but now (pk{x) — x"/^e ^^"^Pkix) where the pk are the Laguerre 
polynomials L^, normalized. 

These processes have scaling limits. If wc scale the Hermite process at the edge we obtain 
the Airy process with corresponding extended Airy kernel [7, 13] 



The Airy process consists of infinitely many curves and as before det (I—K) is the probability 
that no curve passes through at time r^. In the case of greatest interest — {Ck, 
and then the determinant is equal to the probability 



where A{t) is the top curve of the Airy process. (This is what has been called the Airy 
process in the literature. It is convenient for us to use the different terminology.) 

If we scale the Hermite process in the bulk we obtain the sine process with the associated 

extended sine kernel 



If we scale the Laguerre process at the bottom (the "hard edge" ) we obtain the Bessel 
process and its associated extended Bessel kernel 



where ^aiz) = y/zJa{z)^ 

The Airy process A{t) was introduced by Prahofer and Spohn [13] as the limiting sta- 
tionary process for a polynuclear growth model. (See also [7].) It is conjectured that it 
is in fact the hmiting process for a wide class of random growth models. Thus it is more 
significant than the Hermite process. It might be expected that likewise the sine process 
(possibly) and the Bessel process (more likely) will prove to be more significant than the 
unsealed processes. 

For m = 1 the extended Airy kernel reduces to the Airy kernel and it is known [14] that 
then (1.7) is expressible in terms of a solution to Painleve II. It was thus natural for the 

'*Hints of this kernel for m = 2 appear in [10]. 




(1.6) 



Pr {A{n) <Ci,...,A{t^) <U), 



(1.7) 
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authors of [7, 13] to conjecture that the m-dimensional distribution functions (1.7) are also 
expressible in terms of a solution to a system of differential equations. This conjecture was 
established in two different forms, by the authors in [17] and by Adler and van Moerbeke for 
m = 2 in [2].^ Specifically, in [17] we found a system of ordinary differential equations with 
independent variable ^ whose solution determined the probabilities 

Pr {Ar, <Cl + C,---,Ar <U + 0- 

The appeared as parameters in the equations. 

Our first result is a generalization and strengthening of this. We assume that each 
is a finite union of intervals rather than a single interval, and find a total system of partial 
differential equations, with the end-points of the intervals of the Xk as independent variables, 
whose solution determines det {I — K). (When X^ — (^fe, oo) it it easy to recover the system 
of ODEs found in [17].) 

Then we find the analogous systems for the Hermite process (which is more complicated) 
and also for the sine process. Finally we find a system of PDEs for the Bessel process, which 
was the most difficult. It is possible that we could find a system for the Laguerre process also, 
but it would be even more complicated (since Laguerre:Bessel::Hermite:Airy) and probably 
of less interest. 

All of these equations in a sense generalize those for the Hermite, Airy, sine and Laguerre 
kernels foimd in [15], which are the cases when m = 1. Although some of the ingredients 
are the same the equations derived here when m = 1 are not the same as those of [15]. 
For example, the special case of the extended Airy equations for a semi-infinite interval and 
m = 1 is the Painleve II equation whereas in [15] one had to do a little work to get to 
Painleve II from the equations. 

We begin with Section II, where we revisit the class of probabilities for which the correla- 
tion functions were derived in [4] and give a direct derivation of the corresponding Fredholm 
determinant representations. The method has similarities to that of [4] (in fact we adopt 
much of their notation) and the results are equivalent. But we avoid some awkward combi- 
natorics. Our derivation is analogous to that of [16] for random matrix ensembles whereas 
the method of [4] is more like that in [11]. 

In Section III we use the previous result to derive the extended Hermite kernel. This is 
of course not new. But since the derivation does not seem to have seen print before, this 
seems a reasonable place to present it. 

In the following sections we derive the systems of PDEs for the extended Airy, Hermite 
and sine kernels. Presumably the other two could be obtained by scaling the equations for 
Hermite, but Airy is simpler and so we do it first. Moreover all the systems will have the 
same general form, and doing Airy first will simplify the other derivations. 

In Section VII we derive the extended Laguerre kernel, and in Section VIII establish the 

^In [1] the authors had already considered the Hermite process in the case m = 2, in our terminology, 
and found a PDE in r = r2 — ti and the end-points of Xi and X2 for the probability that at time no 
curve passes through Xi. In [2] they deduced for the Airy process by a limiting argument a PDE in ^1, ^2 
and T = Ti — T2 when Xi = (^j, 00). These equations and those we find appear to be unrelated. 
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system of PDEs for the extended Bessel kernel. 



II. Extended kernels 



For the most part we shall follow the notation in [4]. We assume the probability density 
for the eigenvalues Xki {i = I, ■ ■ ■ ,n, k = 1, . . . ,m) is given up to a normalization constant 

by 

P(Aii, . . . , Ai^; . . . ; A^^i, . . . , A^^) 

71 

= n e n detK(Afe,,, Xk+i,j)) A(Ai) A(A^), (2.1) 

k=l k=l 

where Vk and are given functions satisfying some general conditions and A denotes Van- 
dermonde determinant. (Indices i,j in the determinants run from 1 to n, and here Ai resp. 
Xm denotes Xn resp. Xmi-) What we are interested in is the expected value of 

m n 
k=l i=l 

so we integrate this times P over all the Xki- 
We apply the general identity 

j ■■■ j det{ipj{xk))lk=i ■ det{'iljj{xk))lk=i d^^xi) ■ ■ ■ d^^Xn) 



n\ det 



to the integral over An, . . . , Ai„, with the part of the integrand containing these variables. 
This includes two determinants, A(Ai) and the factor det (Hi(Aij, X2j))- The result is that 
this n-tuple integral is replaced by the determinant 



det 



{^j Ale-^(^^) iii(Ai, A2,) (1 + /i(Ai)) dX^ . 



Then we use the same identity to rewrite the integral with respect to the X2i using this 
determinant and the factor det {u2{X2i, X^j)). And so on. At the end we use the determinant 
coming from the previous use of the identity and A(Am)- The end result is that the expected 
value in question is a constant times the determinant of the matrix with i,j entry 



■■■ X\Xle ^=1 llMh:K+i)ll(^ + fk{h))dX,---dX^. (2.2) 

k=l k=l 



By changing the normalization factor we may replace A*^ by any sequence of polynomials, 
which we call Pii(Ai), and replace A^ by any sequence of polynomials, which we call Qmj{Xm)- 

^In our applications fj, will be minus the characteristic function of Xk, so the expected value will equal 
the probability that Xki ^ X^ for all k and i. 
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We choose them so that after these replacements the integral with all the fk set equal to 
zero equals S^. In particular the normalization constant is now equal to 1. 

If we write 

m 

e [[ Uk{Xk: Afe+i) = -E'l2(Ai, A2) -E'23(A2, A3) • • • Ejn-i^rni^m-l: ^m) (2-3) 

fe=l 

(there is some choice in the factors on the right), we see that the matrix in question equals 
the identity matrix plus the matrix with i,j entry 



m— 1 

Qruji^m) dXi ■ ■ ■ dXm- 



/„ m—l 
•■■ I Ph(Ai) Ek,k^i{\kiXk+i) 

k=l 



Ui^ + fkiXk)) - 1 

.k=l 



The bracketed expression may be written as a sum of products, 

XI X /fei(Ai) • • •/fe^(AfeJ. 

r>l k\<--<kr 

Correspondingly the integral is a sum of integrals. Consider the integral corresponding to 
the above-displayed summand. For k > j we define 

Ejk{Xj, Xk) — Ejj^i{Xj, Aj+i) * • • • * Ek-i^k{Xk-i, Afe), 

where the asterisk denotes kernel composition, and set 

Pki{Xk) — j -pLi(Ai) -E'ifc(Ai, Xk) dXi, Qkj{Xk) — j Ekm{Xk, Am) Qmj{Xm) dXm- 

By integrating first with respect to the A^ with k ^ ki, . . . , kr, we see that the corresponding 
integral is equal to 

j " ' j /fci(Afei) Pki,i{Xki) -£'fei,fe2(Afei, Afe2)/fc2(Afc2) • • • 

• • • Ek^_^,kr{Xk^_^i Xk^) fkr{Xkr) Qkr,j{Xkr) dXk^ ' ' ' dXk^- 

We deliberately distributed the / factors as we did since if wc let A^j: be the operator with 
kernel Afe^(Afc, A^) = Eke{Xk, X()f{X() then the above may be written as the single integral 

j fkl{X) Pk^,i{X) Ak-^,k2 ■ ■ ■ ^kr-l,krQkr,j{X) dX. 

(If r = 1 we interpret the operator product to be the identity.) Replacing the index ki by k 
and changing notation, we see that the sum of all of these equals 

/E/.(A)A,.,.)fEE-«,^.,.--.-,.«..(A)lA 

k \r>Oki,...,kr I 

where the inner sum runs over all /c^ > • • • > /ci > A;. (If r = the inner sum is interpreted 
to be Qfe,i(A).) 
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We think of /fe(A) Pk,i{X) as the kth entry of a row matrix and the inner sum 

r>0 ki,...,kr 

as the kth entry of a column matrix. The integrand is the product of these matrices. If we 
use the general fact that det (/ + ST) = det (/ + TS) we see that the determinant of / plus 
the matrix with the above i,j entry is equal to the determinant of / plus the operator with 
matrix kernel having k, £ entry 

n-l / \ 

X XX] ^k,kl^kl,k2 ■ ■ ■ Akr-ukrQkrA^) "Pj? J (/^) /^(/^) , 

3=0 \r>Ok'i_,...,kr I 

where in the inner sum k^ > ■ ■ ■ > k\> k. 

This is the /c, i entry of a certain operator matrix acting from the left on the matrix with 
A;, I entry 

x:q.j(a)p,,,(/.)/,(/.). 

i=o 

That matrix is upper-triangular, all diagonal entries are /, and for /c < £ the /c, I entry equals 

X] ^fe,fel^fel,ik2 ■ ■ ■ ■^kr.l- 

k<ki<—<kr<t 

Elementary algebra shows (even for non-commuting variables A^e) that this is the inverse of 
the upper-triangular matrix with diagonal entries / and k, i entry —A^^ otherwise. 

If we recall that A^^lX, /i) — Ei^^{\ /i) fi{fJ,) then we see that we have shown the following: 
Let H{X, fj,) be the matrix kernel given by 

n-l 
j=0 

let E be the matrix kernel with k,i entry Eke{X, /J.) (thought of as when k > I), and let 
f{li) — diag {fk{iJ>))- Then the expected value equals the determinant of 

7 + (7 - Efy'Hf = (7 - Ef)-' [I+{H- E)f]. 

The factor on the left equals 7 plus a strictly upper-triangular matrix, so its determinant 
equals one. Therefore the expected value equals 

det [7 + (77 -E)/], 

and H — E is the extended kernel. 
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III. The extended Hermite kernel 



We have times ri < ■ ■ ■ < and we set = e^'=^^'=+i, with the conventions Tq = 
—oo, Tm+i = +00 so that Qo = Qm = 0. For the Hermite process the probabihty density is 
given by (1.2) so we are in the case where 

and we want to compute the kernel H — E ol Section II. 
We define the Mehler kernel 

K{q- A,/.) = (7r(l -g^))-V2e-T^^^-T^'^^+T^^^ 
which has the representation 

oo 

K{q- \^x)^Y.^'p,{\)p,{^x)e-^\ (3.1) 

SO 

lK{q;X,iJ,)pi{iJ,)dfx = q'pi{X). (3.2) 

Here Pi are the normalized Hermite polynomials. 

We can write the exponent on the left side of (2.3) as 

m— 1 2 m— 1 -1 

E 'ik \2 _ \ ^ ^ \2 

1 „2^A: 1 _ „2'^fc+l' 

fc=l ^ k=0 ^ 

and so, aside from a normalization constant, the left side of (2.3) is equal to 

m—l 

e"^? n ^i^k] Afe,Ajk+i). 
fe=i 

Thus we may take in (2.3) 

^i2(Ai,A2) = e-^?X(gi; Ai,A2), 
Ek,k+i{^k, ^k+i) — K{qk] Afc, Afc+i), {k > 1). 
It follows from (3.1) that K{q) * K{q') = K{qq') when g, g' > and so 

Eik{X, li) = e"^' K{qi ■ ■ ■ qk-i] A, //). 
In particular we deduce from (3.2) that 

j j pi{X) Eim{\ A«) Pj(a') d^JidX =(?!••• qm-iY j PiW e"^' Pj{X) dA = (gi • • • qm-iYS, 
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Hence we may take 

Pli = Pi: Qmj = (?!••• ?m-l)~^ Pj 

as the polynomials in the previous discussion. We see that 



Pkiil^) = j Pi{X) e K{qi- ■ -qk-i, \li)dX= {qi- ■ -qk-iY Pi{li)e ^\ {k > I), 
QkjW = J K{qk ■ ■ ■ qm-i; A, Ai) Qmj{lj)dix = (gi • • • qk-i)~^ PjW, {k > 1), 



It follows that H is the matrix with /c, I entry 



j=0 



Qi ■ ■■qk-i 



left-multiplied by the matrix diag (e^^ 1 ■ ■ ■ 1) and right-multiplied by the matrix 
diag(l ■■■ ). Similarly E is the strictly upper-triangular matrix with A;,£ entry 

K{qk ■ ■ ■ qe-i', A, fi) left-multiplied by the matrix diag (e^'^ 1 ■ ■ ■ 1). 

Thus we have computed H — E. The actual extended Hcrmite kernel will be a mod- 
ification of this. The determinant is unchanged if we multiply H — E on the left by 
diag(e-^^/^ e'^^/"^ ■■■ e~^^/^) and on the right by diag (e~^^/^ e^^/^ ••• e^^/^). Recalling 
that (pi are the harmonic oscillator functions and recalling the definition of the qk in terms 
of the Tk we see that the expected value in question is equal to 

det[7 + (^-^)/], 

where 

n-l 

and E is the strictly upper-triangular matrix with k, i entry 
If we observe that by (3.1) 

oo 

g(M2_^2)/2 J^(^^r,-r, .X,l_l) = Y^ ^^'-^'^ ipjil^) 

3=0 

when k < i we see that H — E has k, i entry 



{H - E),e = 



"e e^'^^'="'"^V^-(A)(^,(At) ifA;>^, 

j=0 



E e^("'=-"^V,(A)9'i(/^) iik<e, 



j=n 



which is the extended Hermite kernel (1.3). 
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IV. PDEs for the extended Airy kernel 



We consider first the case — (^fe, oo), so that 

det {I-K) = Pr (A(ri) < 6, • • • , ^irm) < U) ■ 

The derivation is simplest here but it will also give the main ideas for all the derivations. 
Observe first that 

dkK^-L5k, (4.1) 

where Sk denotes multiplication by the diagonal matrix with all entries zero except for the 
kth, which equals S{y — ^k)- It follows that if we let R — K {I — K)~^, then 

dk log det(7 - X) = -IV (7 - K)-^ dkK = R^ki^k, ^k)- 

The matrix entries on the right will be among the unknowns. To explain the others, let A{x) 
denote the m x m diagonal matrix diag (Ai(a;)) and X{x) the diagonal matrix diag {Xk{x)), 
where Xk — X(5j.,oo)- Then we define the matrix functions Q{x) and Q{x) by 

Q^{I-K)-'A, Q^Ax{I-K)-' 

(where for Q the operators act on the right). These and R{x,y) are functions of the ^k as 
well as X and y. We define the matrix functions q and r of the only by 

Our unknown functions will be these and the matrix functions q' and q' defined by 

We shall also write and for the matrices {Rxij{Ci,Cj)) i^yijiCiTCj))- 

The ^k are the independent variables in our equations. We denote by the matrix 
diag (^fc) and by the matrix of differentials diag {d^k)- With these notations our system 
of equations is 



dr = -rd^r + d^r^ + ryd^, (4.2) 

dq = d^q'-rd^q, (4.3) 

dq = q'di-qd^r, (4.4) 

H = dCCq-irxdC + dCry)q + dCrq', (4.5) 

dq' = q^d^-q{d^ry + r,dO + q'rd^. (4.6) 



One sees that the right sides involve the diagonal entries of + Vy and the off-diagonal 
entries of and Vy. Wc shall show below that these arc "known" in the sense that they are 

^At points of discontinuity we always take limits from the right. For example we interpret Rij{x, ^j) as 
the limit Rij{x, 
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expressible algebraically in terms of our unknown functions, so the above is a closed system 
of PDEs. 

We begin by establishing the assertions about Tx and Vy. 

In the following D — d/dx, we set p — {I — K)~^ and S — J2k ^k, and r is the diagonal 
matrix diag(rfc). We denote by © the matrix with all entries equal to one. For clarity we 
sometimes write the kernel of an operator in place of the operator itself. 

Lemma 1. We have the commutator relation 

[D, R] = -Q{x)QQ{y) + Rdp+ [r, R\. (4.7) 
Proof. Integrating by parts in (1.6) gives 

[D, K]ij = -Ai(x) ki{y) Xj{y) + Lij{x, Q 6{y - Q + (r^ - r,) Kij{x, y). 
Equivalently, 

[D, K] = -A{x)QA{y) X{y) + + [r, X]. 

To obtain [D, i?] we replace K hy K — I m. the commutators and left- and right-multiply 
by p. The result is (4.7).^ 

If we take the j entry of both sides of (4.7) and set x = ^j, y — we obtain 

i^x + i^y — —qQq + r^ + [t, r] . (4.8) 

Thus all entries of rx + Vy are known. 

For the off-diagonal entries of Tx and r-y we need a second commutator identity. Here M 
is multiplication by x. 

Lemma 2. We have 

[L>'-M, p]^R5px-Ry5p, 
where Ry{x,y) is interpreted as not containing a delta- function summand. 

Proof. We use the facts that — M commutes with L and that M commutes with X. 
These give 

[D^ -M, K]^ [D^ - M, LX] = L [L>2 - M, X] = L [D^, X]^L(5D + D5). 
Using the commutator identity 

[T, (/ - K)-'] = (/ - K)-' [T, K] (/ - K)-\ 
valid for any operators T and K, we deduce 

[D^-M, p]^ pLSDp + pLDSp. 

^Because of the fact pLx = R and our interpretation of Rij{x, ^j) as Rij{x, we are able to write 
R6 p in place oi pL6 p. 
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The first term on the right equals R5 px- The second term equals —Ry5p where Ry is 
interpreted as not containing the delta-function summand. This establishes the lemma. 

Lemma 1 says 

Rx + Ry^ -Q{x) e Q{y) + R5p + [r, R], 
and applying dx — dy to both sides gives 

Rxx - Ryy = -Q'{x) QQ{y) + Q{x) 6 Q'{y) + RxSp - RSpy + [r, Rx - Ry]. 

Lemma 2 says 

Rxx - Ryy - {x -y)R^ R5px - Ry5p. 
Equating the two expressions for Rxx — Ryy gives 

{x - y) R{x, y) = -Q\x) 9 Q{y) + Q{x) 9 Q\y) + {Rx + Ry)6p - RS{px + Py) + [r, Rx - Ry]. 

(4.9) 

Taking the i,j entries and setting x — ^i, y — Cj give 

[^,r] + rrx-ryr^-qQq + qQq+rxr-rry+[T,rx- Ty], 

or 

[r, rx -ry]^q'eq-qeq' + [r, + ry] + r]. (4.10) 

The left side has i,j entry (tj — tj) {vxij — Vyij) and the right side is known.^ Therefore the 
off-diagonal entries of Tx — Vy are known, and therefore so also are the off-diagonal entries of 
Tx and Ty individually. 

To be more explicit we define matrices U and V by 

U = — g 9 g + + [r, r] 

and 

_ {q'Qq-qO q')ij + [r, -qQq+ [r, r] - ^]ij 

n - Tj 

when i ^ j. Then (4.8) says 
and (4.10) gives 



when i ^ j. It follows that for such i,j we have 



dii Tyij d^j Txij = ]j{dii dij)Uij - ]^{dii - dQVij. 



^Here Vxij is notational shorthand for {rx)ij and Vyij for {ry)i 
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The same hold when i — j ii we interpret the second terms to be zero then. More succinctly, 
d^r, + ryd^^ ^{dt U} + ^[dt V], d^Vy + r,d^^ ^{d^ U} - ^[d^ V], 

where the curly brackets indicate anticommutator. These give the explicit representations 
for the terms involving and Vy in the equations. 

With our assertions concerning and Ty established we proceed to derive the equations. 
It follows from the general identity 

dk (/ - K)-' = (/ - K)-' dkK {I - K)-\ 
relation (4.1) and the remark in footnote 7 that 

dkP^-R5kp. (4.11) 
Prom this we obtain (since dkR — dk p) 

dk r ij = dk {Riji^i, Cj)) = {dk Rij){^i, 0) + ^xiji^i, ^j) ^ik + Ryiji^i, 0) ^jk 

^ik ^kj ~l~ Rxiji^ii Cj) ^ik ~l~ Ryiji.^ii Cj) ^jk- 

Multipliying by d^k and summing over k give (4.2). 
Using (4.11) applied to A we obtain 

dk Qij = Q'iji^i) ^ik - {R Sk Q)ij{^i) = Q'iji^i) kk - Tik Qkj ■ (4.12) 

Now multiplying by d^k a-nd summing over k give (4.3). 

It follows from (4.11) that dk Px = —Rx ^k P- Applying this to A gives dk Q' = —Rx ^k Q, 
whose i,j entry evaluated at x — equals —Vxik Qkj- Hence 

dk q'ij = dk Qij (Ci ) = -Txik Qkj + Sik Qij i^i) . (4.13) 

Now we use Lemma 2 again. Applying both sides to A and using the fact that {D'^—M)A — 
we obtain 

Q"{x) - xQ{x) = RSQ' - RySQ. (4.14) 
Taking the i, j entry and evaluating at a; = gives 

Qiji^i) - Qij = - Ty q)ij. 
Substituting this into (4.13) we obtain 

dk q'ij = -Txik Qkj + kk [6 Qij + {fq' - Ty q)ij\. 

Multipliying by d^^ and summing over k give (4.5). 

To obtain the other equations, we point out that identities such as these occur in dual 
pairs. Observe that the function Xj{y) Pjk{y,x) is equal to Xk{x) times pkj{x,y), where p 
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is the resolvent kernel for the matrix kernel with i,j entry Lji{x,y) Xj{y)- Hence Qjk{x) is 
equal to Xfe(x) times the Qkj{x) associated with Lji. The upshot is that for any formula 
involving g or g there is another. We replace q by g* and q with g*. (If a formula involves 
r we replace it by r* and subscripts x and y appearing in r are interchanged.) In this way 
equations (4.4) and (4.6) are consequences of (4.3) and (4.5). 

Let us derive the system of equations found in [17]. We introduce the differential operator 
V = J2k ^k- The system of equations is 

V'^q = ^q + 2qQqq-2[T, r]q, (4.15) 
V'^q = g^ + 2ggeg-2g[T, r], (4.16) 
Vr = -geg + [T,r]. (4.17) 

This can in fact be thought of as a system of ODEs since if we replace ^i, - ■ ■ by 
+ ■ ■ ■ ) Cm + ^ then V — d/d^ and the are parameters in the equations. 

Equation (4.17) follows upon summing over k the coefficients of the d^k in (4.2) and using 
(4.8). Similarly (4.3) gives Vq — q' — rq, so 

"D^ g = P g' + (g g — [r, r]) g — r (g' — r g). (4-18) 

Finally, (4.5) gives 

Vq = -{r^ + Ty) q + ^ q + r q' . 
Substituting this into (4.18) and using (4.8) again give (4.15). We derive (4.16) similarly. 

When m = 1 (4.15) is the Painleve II equation g" = ^g + 2 g^. 

We now consider the more general case where each Xk is a finite union of intervals, 

— i^kl, ^k2) U (Cfe3, Cfe4) U • • • . 

We write dkw for d/d^kw We have 

dky,K= i-ir L6Uy), (4.19) 

where 5kw{y) is the m x m diagonal matrix all of whose entries are except for the kth, 
which equals 6{y — ^kw)- It follows that 

dky, log det(/ -K)^ -Tr (/ - K)-' dk^K = {-lT+^Rkk{^kw, Ckw)- 

The various ^kw are the independent variables. (We shall systematically use u, v and w 
as indices to order the end-points of the intervals of Xj, Xj and X^, respectively.) We now 
define the matrix functions r, g, g, g' and g' of the ^kw by 

and 

^°At points of discontinuity we always take limits from inside Xk- 



15 



These will be the unknown functions in our PDEs. We also define Tx and Vy by 



r. 



x,iu,j 

Observe that r, Vx and Vy are square matrices with rows and columns indexed by the 
end-points kw of the Xk while q, q', q and q' are rectangular matrices. Further notation is 

e = diag(a^), de = diag((-ir+^deto), de = diag(defej, 5 = ^(-ir+^4«,. (4.20) 

k,w 

These are all square matrices but ^, and d$, are indexed by the end-points of the while 

S is m X m. 

With these notations our system of equations is 

dr = -rd^r + d^rx + ryd^, (4.21) 

dq = d^q'-rd^q, (4.22) 

dq = q'd^-qd^r, (4.23) 

dg' = dee?-(r.rfe + ^^er,)g + rfer?', (4.24) 

dq' = ge^e-g(^^er, + r,dO + ^'^^^e- (4.25) 

As before the right sides involve the diagonal entries of Vx+Ty and the off-diagonal entries 

of Tx and Ty, and we must show that these are known. 

It is easy to see that Lemmas 1 and 2 still hold with the new definition of S. Lemma 1 
gives 

'^x,iu,jv ~l~ '^y,iu,jv — ^ ^ Qiu,k 0.1, jv ~l~ ^ ^ ( 1) '^iu,kw '^kw,jv ~l~ ("'"i ''j) fiu,jv 

k,£ k,w 

In matrix terms, 

'"x + '"y = — g © g — r s r -|- [r, r] , 
where s — diag ((— 1)*"+^). Thus Tx + is known. 

What remains is to show that Vx^iujv a-nd Ty^iujv a-re known when iu ^ jv. 
Prom (4.9) we have, using (4.7) again, 

{x - y) R{x, y) = -Q\x) Q Q{y) + Q{x) Q Q\y) + {Rx + Ry)5p - RS{px + Py) + [r, Rx - R^l 

so 

[M, R] = -Q'{x) e (5(y) + Q{x) © g'(y) - Q{x) © g5p(y) + i?5g(x) © Q{y) 

+[t, i?] 5p - R5 [r, i?] + [r, i?^ - i?^/]. 

It follows, as before, that rx,iu,jv ~ ^y,iu,jv is known when i ^ j and so also are rx,iu,jv and 
fy,iu,3v individually. It remains to determine these when i — j but v. 

To do this we use the identity [DM, R]^D [M, R] + [D, R]M to compute 
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[DM, R] = -Q"ix) e Qiy) + Q'{x) O Q'{y) - Q'{x) O Q5p{y) + RJQ{x) O Q{y) 



+[t, i?^] 5p - RJ [r, i?] + [r, i?,;^ - i?,;^,] - yQ{x) 9 g(7/) + yR6p + |/[r, i?]. 

Next we use (4.14), which is the same here. This gives an expression for Q"{x) which we 
substitute into the first term above to obtain 

[DM, R]^-(x + y) Q(x) e Q(y) + Q'{x) O Q'(y) - Q'{x) 6 Q5p(y) + (R, + Ry)5Q(x) O Q{y) 

+R 6 Q'{x) eQ{y) + [r, R^] 6p - RJ [r, R] + [t, R^^ - R^y] + yR6p + y[T, R] . 

The left side equals xR^ + yRy + R and its i, i entry evaluated at {C,iu, C,iv) equals 
Ciurx,iu,iv + ^ivry,iu,iv + riu,iv If we can compute this sum then we know r^^i^^i^ and r^^j^^j^ 
individually since we know rx,iu,iv + ^y,iu,iv and ^ ^iy. To see that the corresponding 
right side is computable observe that the term arising from R^ + Ry is known because of 
Lemma 1, and the diagonal entries of [r, Rxx — Rxy] are zero. Everything else is easily seen 
to be computable except possibly the terms arising from the sum [r, Rx] Sp — RxS [r, R]. Its 
i, i entry equals two times 

Tj y^X~^)^ Rxikjx , Ckw) RkiiCkw, u) — '^{ — ^Y'RxikiXi Ckw) Tk Rki{Ckw , u)- 
k,w k,w 

The two summands corresponding to A; = z cancel. The remaining terms evaluated at 
iCiu, Civ) involve rkw,iv and rx,iu,kw with k ^ i, all of which are known. 

This completes the demonstration that all terms on the right sides of our equations are 
known. This was the hard part. With (4.11) replaced by 

dkwP^ RSkwP, 
the derivation of the equations proceeds exactly as before, and need not be repeated. 

Remark 1. One might wonder whether the systems of equations (4.2)-(4.6) and (4.21)- 

(4.25) are integrable in the sense that one can derive from the equations themselves that the 
differentials of the right sides are zero. Because of the complicated expressions for and 
Ty we have not attempted to show this in general. For equations (4.2)-(4.6), where we have 
relatively simple expressions for the right sides, we verified that this is so when m = 2 or 3. 

Remark 2. We point out how little the equations depend on the operator L, as long as we 
still define K = Lx with X = diag(Xxfc)- Equation (4.21) holds for any integral operator 
L. So does (4.22) if q is defined as before in terms of Q = (/ — K)~^(p, where tp can be 
any function whatsoever. Similarly for q and (4.23). Similarly also for the right hand sides 
of (4.24) and (4.25) except for the first terms d^^q and q^d^. What does depend on the 
specifics of L are the following: 

(i) The expressions for and Ty in terms of the unknowns. We do not see these explicitly 
in the equations. This is where the choice of (p arises. 

(ii) The first terms on the right sides of (4.24) and (4.25), which arise from the computation 
of g". (See (4.14).) 
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All our systems will have the same form as these, most of the equations being universal, i.e., 
independent of the specific L or ip}'^ In most cases there will be two functions such as </?. 
That will add to the number of equations but not their complexity. The main difficulty in 
all cases will be (i). 

V. PDEs for the extended Hermite kernel 

We modify (1.3) by setting 



Lij{x,y) 



e(^-") (-^--^) ^x) ifkiy) if ^ > J, 

A;=0 



_ ^ e^^-^^^-^--^^^k{x)My) if ^<J. 

k=n 

The extra factors e"'*^'^*"'^-'^ do not change the determinant. 

Again we consider first the case where X^. = ({fe, oo). We define R and p as before, and 
again 

9felog det(/ - K)^ Rkki^k, Cfe)- 
Now we shall have more unknown functions. We set 

99 = (2n)i/Vn, ^ = (2n)i/Vn-i, 

and define 

Q = pip, P = pip, Q = ifXp, P = ipXp. 

Our unknowns will be, in addition to r^j = Rij{Ci, ^j), the matrix functions q, q, p and p 
given by 

and 

olj-Qiji^i): Qlj-QiAQ: i^ii-Pim^ i^ij-HM- 

Again ^ denotes the matrix diag (^fc) and denotes diag {d^k)- 
With these notations our system of equations is 



dr = 


-rd^ 


r + d^r^ + Ty d^, 


(5.1) 


dq = 


diq'- 


-rd^q, 


(5.2) 


dq = 


q'd^- 


- qd^r, 


(5.3) 


dq' = 


dae 


-2n- l)q- (fx d^ + d^ Ty) q + d^r g'. 


(5.4) 


dq' = 


ne- 


-2n-l)di-q{diry + rx di) + q' r d^. 


(5.5) 


dp = 


dip - 


-rd^p, 


(5.6) 


dp = 


p'di- 


-pdi r. 


(5.7) 


dp' = 


dae 


— 2n + 1) p — ij-x d^ + d^ Ty) p + dC, rp', 


(5.8) 


dp' = 


Pie- 


- 2n + 1) d^ ~ p [d^Ty + Tx dC,) + p r d^. 


(5.9) 



'^This splitting into universal and nonuniversal equations was also a feature of [15]. 
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By Remark 2 and duality (each equation for q or p giving rise to one for q or p) all we 
have to show is that the diagonal entries of + Ty and the off- diagonal entries of and Vy 
are known (i.e., expressible in terms of the unknowns) and to derive (5.4) and (5.8). 

We begin by finding a substitute for Lemma 1. We write for D ± M. 
Lemma 3. We have 

D+Lij - e^^-^^ UjD+ = -^P{x) ip{y), e^'"^^ D-L^j - L^jD- = -ip{x) ^P{y). (5.10) 
Proof. Let J be the operator on i>^(R) with kernel 

n-l 
k=0 



and set ak = yk/2. We have the formulas 

X(pk = Ofc+l (Pk+i + dk Ifik-l, f'k = -dk+i <fk+l + (Pk-i- 



Therefore 



n-l 



{x + dx) J{x,y) = 2 ^a^ak^k-i {x) (fik iv) , 

k=0 

{y - dy) J{x,y)^2 ^ a'' ak+i <Pk{x) <Pk+i{y)- 



k=0 



This gives 

[{x + dx) - a {y - dy) ] J{x,y) = -2 a"" Un (pn-iix) <fn{y)- 

If we take a = e^^^^^ and multiply by e^"*-'^*"'^^'' we obtain the first identity of (5.10) when 
i > j. If 0" < 1 and one takes n — > oo in the last identity for J one gets zero for the right 
sides. It follows that replacing Ylk=o ~ Sfcln definition does not change the right 

side. Thus wc obtain the identity for i < j as well. The second identity of (5.10) is obtained 
from the first by taking adjoints and using the fact that Lij is self-adjoint. 

We can now find the analogue (actually, analogues) of Lemma 1. Observe that since 
r = diag (r^) we have e'^ — diag (e"^'). 

Lemma 4. We have 

e-^D+ R - Re-^D+ = -P{x) e"^ 6 Q{y) + i?5e->, (5.11) 

e'D- R - Re^D- = -Q{x) 6 P{y) + Rde" p, (5.12) 
Proof. If we multiply the relations (5.10) on the right by X and use the fact [D'^, x\ — 5 we 
obtain 

e-^D+ K-K e-^D+ = -e'^ ij;{x) Q x{y) ^{y) + Lde'^ , 
e'^D- K-Ke''D- ^ -e^ ^p{x) O X{y) V'(y) + Lde" . 
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We replace K on the left hy K — I and left- and right-multiply by p, and the result follows. 
(We used the fact that e^"^ commutes with the matrix functions </? and '0.) 

If we take i, j entries in (5.11) and (5.12) and set x = ^j, y = we obtain 

e^^r^ + Ty er'^ = — e^^ + r e^^ ^ — pe~'^ Qq + r e~'^ r, 

Tx + Ty = — r — qe^ Qp + r r. (5.13) 

The right sides here are known. If we add and subtract these identities and take i, j entries 
we obtain 

2 (cosh Ti Txij + cosh Tj Tyij) = ■ ■ ■ , (5.14) 

2 (sinh Ti r^ij + sinh tj Vyij) = ■■■ , (5.15) 

where the dots on the right represent known quantities. The first relation with j = i gives 
r^ii + Tyii. If the two relations are thought of a system of equations for r^ij and r^y the 
determinant of the system is nonzero when i ^ j. Therefore we can solve for r^ij and ryij 
individually then. 

What remains is to derive (5.4) and (5.8). For this we need the analogue of Lemma 2. 
Lemma 5. We have 

[D^-M"^, p]=R5p:,-Ry5p, (5.16) 
where Ry{x,y) is interpreted as not containing a delta- function summand. 

The proof is analogous to that of Lemma 2. Here we use the fact that D"^ — M"^ commutes 
with L, a consequence of the fact that each ipk is an eigenfunction of — M^. 

Since (p is an eigenfunction of — with eigenvalue — 2n — 1 and ip an eigenfunction 
with eigenvalue —2n -\- 1 applying both sides of (5.16) to (p and to ifj gives 

Q" -x^Q + {2n+l)Q = R6Q' - Ry6Q, (5.17) 

P" -x^P+{2n-l)P = R6P' - Ry6P. (5.18) 

We have (4.13) here just as before. Taking the i,j entry in (5.17) and evaluating at 
x^^i gives 

Qiji^i) - i^i - 2n - 1) qij ^{rq' - Ty q)ij. 
Substituting this into (4.13) we obtain 

dk 4j = -Txik Qkj + Sik [{^i - 2n - 1) qij + {rq' - Ty q)i-[ , 
which is (5.4). Equation (5.8) is established in exactly the same way using (5.18). 

We can also derive a system analogous to equations (4.15)-(4.17): 



V^q = {^'^ -2n-l)q-2Vr-q, (5.19) 

V^q = q{^'^ -2n-l) -2 q-Vr, (5.20) 

V^p = {(^ - 2n + l)p-2Vr -p, (5.21) 

V'^p = pi^"^ -2n + l) -2 p-Vr, (5.22) 

Vr = -r'^ + r^ + r-y. (5.23) 
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These equations are not as simple as (4.15)-(4.17) since the expressions for the entries of 
+ Ty are messy. The last equation we already know. The other equations are derived as 
for Airy: Summing the coefficients of d^k in (5.2) gives V q — q' — r q, so 

V^q^Vq' -V r -q-riq' -rq). (5.24) 

Similarly (5.4) gives 

V q' = {^"^ — 2n — 1) q + r q' — (r^, + Vy) q. 

Substituting this into (5.24) and using (5.23) again give (5.19). We derive (5.21) similarly, 
and (5.20) and (5.22) are obtained by duality. 

In case m — 1 (5.13) gives + Vy — r"^ — pq, and our system of equations becomes 

r'^-pq, q"^ie-2n-l)q + 2q^p, p" ^ - 2n + 1) p + 2p' q. 

From the last two we find {pq' — qp')' = pq" — qp" = —2pq, and by the first equation this is 
2r'. Thus pq' — qp' = 2r. Using this, and successively computing r", r'" and r"" using the 
differentiation formulas, we arrive at 

r"" = 4(f - 2n)r" + 4^r' - 12rV" - 4r - (4(^2 - 2n)r')' - 4(^r)' - 6(r' y, 

and so 

r'" = 4(e'-2n)r'-4er-6r'2. 
This is the third-order equation found in [15] which integrates to Painleve IV. 

We turn to the more general case where each is a finite union of intervals, and will 



again use the notations (4.20). The equations are 

dr — —rd^r + d^Tx + Tyd^, (5.25) 

dq = d^q'-rd^q, (5.26) 

dq = q'd^-qd^r, (5.27) 

dq' = di{i^ -2n-l)q- {r,^d^ + d^ry)q + dirq', (5.28) 

dq' = qi^^"^ -2n-l)di-q{diry + r^di) + ^rdi. (5.29) 

dp — d^p' — rd^p, (5.30) 

dp = p'd^ — pd^r, (5.31) 

dp' = di{i'^ -2n+l)p- {r^d^ + d^ry)p + dirp', (5.32) 

dp' = p{^'^ -2n+l)d^-p{d^ry + r:cdCj+p' rd^. (5.33) 



Nothing is new here except to establish that the terms involving Vx and Vy on the right 
are known. As usual those that occur are the diagonal entries of + Vy and the off-diagonal 
entries of and Vy. In our case the terms Vxij and Vyij in (5.14) and (5.15) are replaced 
by Tx^iujv and Vy^iujv and the relations show that these are known when i ^ j and that 
the rx,iu,iv + fy,iu,w are known. It remains to show that rx,iu,w and ry^iu,iv are known when 
u ^ V. 
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Prom (5.12), which says 

[e^D", R] = -Q{x) 9 P{y) + R6e^p, 

we deduce 

[{e^D-f, R] = e^D- {-Q{x) 9 P{y) + RSe" p) + {-Q{x) 9 P(2/) + R5e^ p) e^D'. 

We use S = T ioi matrix functions S and T to denote that the differences Siu^i^iCiu, dv) ~ 
Tiu,iv{^iu, ^iv) are known. If we keep in mind that q, q', p and p' are among our unknowns, 
we see that it follows from the above, after multiplying by e"^"^, that 

[D^ - 2MD, R] = RJe'^pe-^ - RSe'^pyC'^ 

If we subtract this from (5.16) we obtain (since [M^, R] is known) 

2[MD, R] = {R5p^ - RJe'^pe-^) + {RSe" pye'^ - RySp). (5.34) 

Consider the first term on the right. Its iu, iv entry evaluated at {^iu, ^iv) equals 

^ ^ ^w. kw{ 1) fx,kw,iv ^ ] fx^iu, kw ( 1) ''1"kw,iv(^ 
kw kw 

The terms of both sums corresponding to k ^ i are known. So remaining as unknown is the 
sum 

^ 1) (j'iUjiw '^XjiWjiv x,iu,iw'^' iWjiv) ■ 
w 

Analogously the second term on the right of (5.34) is a known quantity plus 

^ X (j'iUjiw ^y,iw,iv ^y,iu,iw ^ iw,iv) ■ 

w 

Adding this to the last sum gives 

But this is known since, as we saw at the beginning, the r^j^.tv + ^y,iu,iv are known. 

We have shown that [MD, i?] = a known matrix function. Its iu, jv entry evaluated at 
{^iu: ^iv) equals 

^iu'^x,iu,jv ^~ ^iv '^y,iu,jv ~l~ '^y,iu,jvi 

^iu^x,iu,jv ~^ ^iv ^y,iu,jv known. But SO is Tx^iu^jv ~^ ^y^iu^jv' Therefore Tx^iu^jv and Ty^iu^^ju 
are both known when u^v. 
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VI. PDEs for the extended sine kernel 



If wc make the substitutions Tj — >• Tijln, x — > x/\/2n, y y/\/2n in the extended 
Hermite kernel (1.3) and let n — > oo we obtain the extended sine kernel 



Lij{x,y) = < 



[-1 



(.n-Tj) 



cos z{x — y)dz if i > j, 



= (Ti-Tj) 



COS z{x — y)dz ii i < j. 



Here we set 



ip{x) = sinx, ipi^x) = cosx, 



and then the other definitions are exactly as in Hermite with the above replacements. The 
unknowns now will be only r, q, q, p and p and the equations for general are 

(6.1) 
(6.2) 
(6.3) 
(6.4) 
(6.5) 

(Recall that s = diag ((-1)"'+^) .) 

We know that equation (6.1) is completely general, as are the equations 



dr — 


—rd^r + d^ 


t + Ty d^, 


dq — 


{p + t^sq) — 


rdiq, 


dq — 


{p + qsr) d^ - 


qd^r, 


dp = 


d^ (— g + rsp) 


-rd^p, 


dp = 


{—q + psr) d^ 


— pd^r. 



dq 


= d^q' 


-rd^q, 


dq 


= q'k 


- qd^r, 


dp 


= kp' 


-rd^p, 


dp 


= p'k 


— pd^r. 



To derive (6.2)-(6.5) from these we establish the formulas 

q' — p + rsq, p' = —q + rsp, q' = p + qsr, p' = —q + psr. 

We have [D, L] = 0, whence [D, K] = L6, whence 

R^ + Ry^[D, p]^RSp. 

Applying (6.7) on the left to (p and ip, using (p' — ip, ip' — —(p, we obtain 

Q'{x)-P{x) = iRdQ){x), P'{x) + Q{x) = {RSP){x). 

Since q^^j — Q\j{^iu) and p ^ = i^j(^i„), the first two relations of (6.6) follow, and the 
others are analogous. 



(6.6) 



(6.7) 
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So equations (6.1)-(6.5) hold, and it remains to deal with the entries of rx and ap- 
pearing on the right side of (6.1). We have to show that the diagonal entries of + Vy are 
known and that rx,iu,jv and Vy^iujv are known when iu ^ jv. 

It follows from (6.7) that + ry — rsr, and so all entries of the sum are known. 

Next, for i > j we have 

2 (Tj - Tj) Lxij = -2 (Tj - Tj) I (""'-""^^ z sin z{x - y) dz 

J 

= — e'^*"'^-'' sin(a; — y) + {x — y) Lij. 
The same holds when i < j. Since [tD, L]ij — ti L^ij + Tj Lyij and Ly — —Lx, this gives 

2 [tD, L] = -e^ ® e"^ sin(a; - y) + [M, L], 

where (g) e~'^ is the matrix with i,j entry e'^'~'^K (This is not a tensor product.) Hence 

2 [tD, K] = -e^ e"^ sin(a; - y) X{y) + [M, /sT] + 2 K5t. 

Replacing K hy K — I m. the commutators and applying p left and right give 

2 [tD, i?] = P{x) ® e"" Q{y) - Q{x) ® e"" P(|/) + [M, i?] + 2 R5Tp. 

The i,j entry of the left side evaluated at (,^j„, S,jv) equals twice TiVx^iujv + 'Tj^y,iu,jv ■, so 
these are known. We deduce, since Tx^iujv + ^y,iu,jv is known, that rx,iu,jv and r-y^i^jv are 
both known when i ^ j. Just as before, the trickier part is to show that rx,iu,iv and Vy^iu^i^, 
are known when u ^ v. 

We compute 

[DM, i?] = [D, i?] M + D [M, i?] 

^R5py + 2 (dx [tD, R] - Rx5tp) - P'{x) ® e'^ Q{y) + Q'{x) ® e'^ P{y). (6.8) 

The i,i entry of [tD, R\ evaluated at (^j„, ^jy) equals Ti{Rx + Ry){im, iiv)- Hence, since 
dx{Rx + Ry) — Rx^P by (6.7), the i,i entry of 9^ [tD, R\ evaluated at (^i^, ^j^,) equals 
Ti{Rx5p)ii{^iu, Civ)- It follows that i,i entry of [ri?, R] — RxSrp evaluated at {^iu, Cjv) 
equals 

^ X -i) (^i "^k) '^x,iu,kw ^kw,iv 
kw 

Since we need only sum over k ^ i all these terms are known. So are the other terms of (6.8) 
evaluated at {^iu, Cjv)- 

We have shown that the i, i entry of [DM, R] evaluated at (^i^, ^jy) is known. This equals 

'^iu.iv ~l~ du '^x,iu,jv ~l~ dv '^y,iu, jv ThuS '^x,iu,jv ~l~ dv '^y,iu,jv is kuOWU. SiuCC Tx^iu,jv ~l~ ^y,iu,jv 

is known and ^ Civ so also are rx,iu,jv and Vy^iujv known. 

Let us see what these give in the case m = 1 for a single interval (—t, t). Here 
^1 = —t, ^2 = t. If we use the fact that K{—x, —y) = K{x, y) and the evenness of cosine and 
the oddncss of sine wc get q2 = —qi, P2 = Pi and if we use the fact that R{x, y) = R{y, x) 
for X, y e {-t, t) we get Vij = Vji. 
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We use the notations r — ru, r — ru- If we observe that d/dt — 82 — 81 then (6.6) gives 



dqi ^ _ dp 

and (6.1) gives 



aqi ^ _ api 



— = r +r -r^c-ry 

and the trivial relation df/dt = —f^ + fy. The general relation + Vy — rsr gives in the 
present notation + Vy = r"^ — f"^, and so 

dr o 

Finally (6.8) gives 

f-tf,+tfy= -Pii-t) Q2{t) + Q',{-t) P2{t) = -f^iQli-t) Pl{-t)). 

Thus 

d , d . . 

^(^r-) = --(g,p,), 

which gives 

QiPi 
r — . 



VII. The Laguerre process 

The Dyson process r — > A{t) on the space of p x n complex matrices (we assume p > n) 
is specified by its finite-dimensional distribution functions. The probability measure on 
Ak — A{Tk) (/c = 1, . . . , m) is a normalization constant times (1.4), which may be written 

X fj exp ( Tr (A*A,_i + he) J dAi--- dA^. (7.1) 

j=2 \^ Qj J 

(Here "he" is an abbreviation for "Hermitian conjugate".) We show how to derive (1.5) from 
this. 

Any complex matrix p x n complex matrix A can, by the singular value decomposition 
(SVD) theorem, be written as 

A = UDV* 

where U is a. p x p unitary matrix, V is an n x n unitary matrix and D is a p x n matrix 
all of whose entries are zero except for the diagonal consisting of the singular values of A. 
Thus we write each Aj as 

A, = u,D,v; 
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with the goal of eventually integrating over the unitaries Uj and Vj. Of course, 



Tr:{A*Aj)^Tv{D*Dj)^J2X, 



n 



where Xj — dijj. 

Let us examine one term 

Tr [A*Aj^i + hc) 

appearing in the exponential of the second product in (7.1). Using the SVD representation 
we have terms 

IV (Vj D] u; v*_, + he) 

The integrals over the unitary group (Haar measure) are both left- and right- invariant. Thus 
in the V^_i integration we let 

V* ^ ^v* , V* 

so that the trace term becomes 

IV {D*U;U,^^D,^^V;_, + hc). 
In the integration we let UjUj-i and the trace becomes 

IV {d* u,., d,_, v;_, + he) . 

Thus, we have integrals of the form 

/ / ^""p (t^ ^-^ + ^'^) '^''^^'-'^ '^''^^'-'^ 

Let S denote an nxp complex matrix, T apxn complex matrix and U (resp. V) elements 
of the unitary group of pxp (resp. n x n) matrices. We assume p > n and set a = p — n. We 
let Si resp. ti denote the eigenvalues of SS* resp. T*T. The Harish-Chandra/Itzykson-Zuber 
integral for rectangular matrices (see, e.g., [18]) is 



J J exp (cTr {SUTV* + he)) di^{U) dniV) 



Cp,n,c det (/„(2c^a,6, 



A(a)A(6) YlUisiUyi^ 



Here c can be any constant, a — p — la is the modified Bessel function and Cp,n,c is a 
known constant. 

When the qj — the measure (7.1) must reduce, after integration over the unitary parts, 
to the well-known Laguerre measure. It follows that (7.1) becomes after integration over the 
unitary parts a normalization constant times 
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(2 \ n 

k=l k=l \ \^ %+l )) 

n n 

X A(Ai) A(A^) n Aif n >^m' ^^11 • • • d^mn, 

i=l 1=1 

which is (1.5). 

Wc shall now compute the extended kernel using the method of Section II. This density 
is not quite of the form (2.1) because of the last factors in the integrand here. Consequently 
in (2.2) there are extra factors Af^^ and A^/^, and so in the discussion that follows Pii(A) 
and Qmi{X) are no longer polynomials of degree i but A"/^ times polynomials of degree i. 

We have now 

We introduce the Hille-Hardy kernel (the analogue of the Mehler kernel) 
which has the representation 

oo 

i^(g;A,/.) = E?'>?WK(/^))""e-^ 

1=0 

where pf are the Laguerre polynomials Lf, normalized. It follows that 

fOO 

/ K{q;X,fi)pnfi)dfi = q''pnX) (7.2) 
Jo 

and so again K{q) * K{q') — K{qq'). 
Now we may take in (2.3) 

Ei2(Ai, A2) = e-^i (Ai/A2)"/'ir(gi; Ai, A2), 

Ek,k+i{^k, ^k+i) — (Afe/Afe+i)"'^^ir(g'fe; Afe,Afe+i), {k > 1), 

and so 

E,k{X, ix) = e-^ (A//x)"/2 K{q^ ■ ■ ■ q,_,; A, //). 
We deduce from (7.2) that 

/ / X"/'pmE,^{\pi)pL"/'p'^{pi)dpidX 
= (gi • • • qm-if^ J A>nA) e-"pJ(A) dX^iqi--- q^-if^ S^j. 
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Hence we may take 



We see that 



Pkiif^) = / A"/2 Pu{X) e-' (A///)"/2 K{q, . . . q,_,; A, //) dX 

= (gi---gfc-if pr(/i)/^"/'e-^ {k>l), 

QkjW = J (A//^)"/' K{qk ■ ■ ■ Qm-i; A, //) g„.,(//) A^"/' d// 
= (?i---?fc-i)-'^A'^/2pJ(A), (A;>1), 



It follows that H is the matrix with k, £ entry 




Pj{X)pj{fj,) 



left-multiphed by the matrix diag {X'^/'^e"^ A"/^ • • • A"/^) and right-multiphed by the matrix 
diag (/i"/^ /i°/^e~^ ••• //"/^e~^). Similarly is the strictly upper-triangular matrix with /c, £ 
entry (A//i)"/^ -f^(gfc • ■ • g^-i; A, fi) Icft-multiphed by the matrix diag (A"/^e~^ A"/^ ■ ■ ■ A"/^). 
Now we use the fact that the determinant is unchanged if we multiply on the left by 
diag (e'^/^ e"*^/^ • • ■ e"^/^) and on the right by diag (e^^/^ e''/^ • • • e^/^). 

In this way we find the analogue of the kernel which was denoted by H—E in Section 111. It 
is now given by (1.3) but with coefficients e^*^^'^'"'^-') and with (pk{x) equal to a;"/^ e~^^'^ p'^{x). 
This is the extended Laguerre kernel. 



If we make the substitutions Tj — > Ti/2n, x — > x^/4n, y — > y'^/A.n in the extended 
Laguerre kernel and then let n — > oo we obtain the extended Bessel kernel 



Let us immediately explain the difficulty. In the previous cases we were able to find one 
commutator for L involving D and another involving D^, the latter arising from the differ- 
ential operator whose eigenfunctions appear in the integrand or summand of the expression 



VIII. PDEs for the extended Bessel kernel 




where 
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for the kernel. (For the extended Airy kernel these were given in Lemmas 1 and 2.) These 
enabled us to express and Vy in terms of the unknown functions. 

Here there does not seem to be a commutator involving the first power of D. We are able 
to find two relations involving the first power of but each involves both a commutator and 
an anticommutator. Fortunately we are able to deduce from these relations three commutator 
relations involving D^, and these relations will enable us to show that the derivatives of 
and Ty are expressible in terms of and Vy and the other unknown functions. The upshot is 
that we are able to find a system of PDEs in which Tx and Vy are now among the unknowns. 
Although the system of equations seems no more complicated than those we have already 
derived (just larger) it is actually much more so because of the expressions for the derivatives 
of Tx and Ty in terms of the unknown functions. 

To state the equations, we define (p and if} by 

From these we define q in the usual way. But now we set 

P ^ {I - Xy^Mi/j, P ^ M^X{I - K)-^, 

and from these we define p and p in the usual way. (The reason we do this is that eventually 
it is these p and p which will arise in the expressions for the derivatives of Tx and Ty.) With 
these notations our system of equations, in the general case where each Xj^ is a finite union 
of intervals, is 



dr = —rd^r + d^Vx + ryd^, (7.1) 

drx = -Tx d^r + d^ r^x + Txy d^, (7.2) 

dry = -rd^ry + d^rxy + Tyyd^, (7.3) 

dq = diq'-rdiq, (7.4) 

dq = qd^-qd^r, (7.5) 

dq' = di{{a^-\)i-\-q)-{rxdi + diry)q + dirq', (7.6) 

d~q' = {{a'-\)q^-'-q)k-q{d^ry + rxdO + q'rd^. (7.7) 

dp = dC^p' — r d^ p, (7.8) 

dp = p d^ — p dC^ r, (7.9) 

dp' = di{{a^ -\)i-^p + 2q-p)-{rxdi + diry)p + dirp', (7.10) 

dp' = {{a^ -\)pi-^ + 2q-p)di-p{diry + rxdi)+p'rdi. (7.11) 



Equations (7.2) and (7.1) are obtained in the same way as (4.21). We have 

9k'w'^x,iu,jv — dkvjRxiji^^iu-i ^jv) — ( 1) Rxiki,Ciuj Ckw^ Rkji,Ckwj Cjv) 
~^Rxxij{,dui Cjv)^iu,kw ~l~ RxyijiCiui Cjv^^jv,kw 

This gives (7.2) and (7.3) is analogous. 
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So all the equations are universal except for (7.6) and (7.10) and their duals. What we 
have to do is show that the diagonal entries of r^x + f^y and r^y + Vyy, and the off-diagonal 
entries of Txx, Txy and Vyy are all known, and to establish equations (7.6) and (7.10). 

To begin, we denote by the kernels where $^(2^^) ^aiyz) in the integrand is replaced 

by 

When a = —1/2, L"*" is essentially the extended sine kernel and some of the formulas we 
derive here will specialize to those obtained in Section VI. We use the notations (5 — j + a 
and 

C{x, y) = (p{x) i){y) - ■4j{x)(p{y), ri{x, y) = (p{x) ij{y) + ^(x) if{y), 

n = (e("^-"^)/2). 

After integration by parts and some computation using the differentiation formulas 

$'Jz) = -$„+i(;^) +/?^-i$«(z), $:,+i(z) = <^>^{z)-(3z-^<^^+^{z) 
we find that 

L+ = C(x, y) + f^L- + -^{x - y) L+, 

Ti Tj X Ti Tj 

= —na^, y) + ^L- + -^{y - x) L+, 

n - Tj y Ti- Tj 

L- = r]{x, y) + ^L+ + —^{x + y)L-, 

Tj Tj X Ti Tj 

Ly = v{x, y) + -L+ + -^—{x + y)L-. 

n-rj y n-Tj 

Here the i,j entries of the matrices L"^ and Q are to be understood. 

If we add the first two identities and subtract the last two we obtain we obtain the 
commutator-anticommutator pair 

[D,L+]=P{M-\L-}, {D,L-} = P[M-\L+]. (7.12) 

To obtain another pair, first multiply the first two identities by Ti — Tj and subtract, 
getting 

in - Tj) (L+ - L+) = 2 Q C(x, y)+{^- in - Tj) L- + 2{x - y) L+. 
Using the first two identities again we can write the left side as 
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^ \ n-Tj y n-Tj 



Thus 



n L+ + T,- L+ - I^^Ti + ^T^j L- + {x-y)L+ + n ((x, y). 
L+ + Tj L+ = + ^Tj^ L- + {x-y)L+ + n ({x, y). 



In other words 

[tD - M, L+] = P{M-W, L-} + nC{x,y). 
Next multiply the last two identities by — Tj and add, getting 

in - Tj) (L- + L-) = -2 Q 77(x, y) + + ^) {n - r,) + 2{x + y) 
The left side may be rewritten 

n L~ - Tj Ly + Ti ( — n rj{x, y) + -L+^ — {x + y) L') 

\ n-Tj y Ti- Tj J 

-Tj ( —nri{x,y) + -L+ + —^{x + y)L-] 



Ti Tj X Ti Tj 



Thus 



Ti - Tj Ly + ^^Tj - j L+ + {x + y)L - Q rj{x, y). 
'i L~ - Tj Ly = (^Ti - ^T^ L+ + (x + y)L- -n r)(x, y). 



In other words 

{tD - M, L-} = /9[M-V, L+] - n7]{x,y). 
Thus we have our second commutator-anticommutator pair 

[tD - M, L+] = /3{M-V, L-} + n Cix, y), 
{tD - M, L-} = (3[M-^T, L+] - n 7]{x, y). 
Now we have the following. 
Lemma. Suppose A and B are such that 

[A, L+] = {B, L-} + F, {A, L-} = [B, L+] + G. 
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Then 

[A' - B\ L+] = [[A Bl L-\ + {A, F} + {S, G}, 
[A' - B\ L-] = [[A, S], L+] + [A, G\ + [S, F]. 

Proof. We have 

[A\ L+] = {A, [A, L+]} = {.1, {S, L-}} + {A, F}. 
By the general identity 

{A, {B, C}} = [[A, S], C] + {B, {A C}} 

the first term on the right side above may be written 

[[A, B],L-] + {B, {A, L-}] = [[A, B], L"] + {B, [B, L+]} + {B, G} 

^[[A,BlL-] + [B\ L+] + {S, G}. 
This estabhshes the first stated identity. 
For the second we write 

[A, L-\ = [A, {A, L-}\ = [A, [B, L+]] + [A, G]. 

By the general identity 

[A, [B, C]] + [B, [C, A]] + [C, [A, B]] = 

the first term on the right side above may be written 

-[B, [L+, A]] - [L+, [A, B]] = [B, {B, L"}] + [B, F] - [L+, [A, B]] 

= L-] + P, S],L+] + [S, F]. 

This gives the second identity. 

We have obtained in (7.12) and (7.13)-(7.14) two quadruples [Ai, Bi, Fi, Gi) and 
{A2, B2, F2, G2) satisfying the hypothesis of the lemma. Each gives commutator relations 
involving L+ and L" . However [A-i + A2, Bi + S2, Fi + F2, Gi + G2) will also satisfy the 
hypothesis of the lemma and so gives commutator relations involving L"*" and L~ . If we 
subtract from these the relations resulting from the other two we obtain 

[AA'+A'A-BB'-B'B, L+] = [[A, B']+[A' , B], L-]+{A, F'}+{A', F}+{S, G"}+{S', G}, 

[AA! + AA-BB'-B'B, L"] = [[A, B'] + [A', B], L+] + [A, G'] + K G] + [B, F'] + [B\ F]. 

So in the end we will obtain three pairs of commutator relations involving L"*" and . If 
we add the identites in each pair and divide by 2 we obtain three commutator identities for 
L. For the explicit computations we have to keep in mind that all matrices and operators 
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commute with (/? and ijj, and D and M commute with r and Q. We write down the results, 
sparing the reader the details: 

[D^ + p{l-p) M-\ L] = 0, 
2 [tD^ - MD + - /9)rM-^ L] 
^Vl(^P®^1){D - I3M-^) - {D + (3M-^)i/j®ip) , 
[(r D-Mf + P{l-P) M'^, L] 
^ n^p ® {t D - M - (3 T M-^) - {tD - M + PTM-^)n'ilj®(p. 

The differentiation formula for ^a+i is in our notation {D+P M~^) ip — (p. Also, an operator 
acts on (p <^ip from the right by applying its transpose to ■0. Using these facts we see that 
the last two identities simplify to 

[tD^ - MD + (3(1 - (3)t M-^, L]^ -2nip0ip, 

[{rD-My + l3{l-f3) T^M-^, L] = -2VLt p ® ^ - ® Mip - Mip ® 

The commutator identities for L lead as before to commutator identities for K — LX- 
They are 

[D^ + M-\ K] = L {6D + D5), 

[tD^ - MD + (3(1 - f^)T M-^ K] = -2Q(p ^ (fX + L (r {6D + DS) - MS) , 
[{t D - Mf + 13 {I - I3)t^ M-^ , K] = -2nT p ® ^X - ® M^X - M^ ® ^X) 

+L (t^ {6D + D6) -2tM6) . 

We are ready to apply p = [I — K)~^ to both sides. The only functions that appear on 
the right sides are </? and M'^, which is why we define 

g = (7-x)-V, Q^<px{i -K)-\ p^{i -Ky^M^ij, p^MiJX{i-Ky\ 

Then we deduce 

[D^ + M-^, R]^ RSp^- RySp, (7.15) 

[tD^ - MD + ^{1-(3)t M-^ R] ^ -2Q(x)nQ(y) + Rr5 p^ - Ryr5 p - RC p, 

[(r D-Mf+P{l- (3y M-^ R] = -2Q{x) Q{y) - Q{x) n P{y) + P{x) nQ{y) 

+R px — Ry t'^S p — Rt^ p. 

We now show that the diagonal entries of r^x + ^xy and Vxy + Vyy, and the ofi^-diagonal 
entries of Vxx, ^xy and Tyy can all be expressed in terms of the unknowns. 

Wc use the symbol = here to mean that the difference of the quantities on its left and 
right is expressible in terms of Q, -P, Q, P and i?, but no derivatives of these functions. 
The three commutator identities above yield in this notation the relations 

Rxx - Ryy = RSpx - RySp, (7-16) 

n Rxx - rj Ryy = xRx + y Ry + Rr6 px - Ry tS p, (7.17) 

T^Rxx-T^Ryy = 2 T X Rx + 2 y Ry T + Rt^6 Px - Ry p. (7.18) 
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Consider first the case i ^ j- It follows from any pair of the above equations (everything 
now is to be evaluated at {^iu,Cjv)) that both R^x and Ryy are known. If we call the right 
sides above A, B and C then 

' I I A 

Ti Ti B 



T-2 r 

Tj U 



= 0. 



If we differentiate with respect to x we deduce that the sum of all terms involving R^y is 
known. (Since our unknowns involved up to one derivative, this is why in our definition of 
= we required that no derivatives were involved in the difference.) This sum is 

-TiTj{Tj - Ti) RxySR - {Tj - Ti) {yRxy - RxyTSR) + {tj - Ti) {2Tjy Rxy - RxyT^SR). 

Dividing this by tj — Ti, evaluating at {^iu, ^jv) and expanding we obtain 



xy,iu,jv 



k.w 



The terms involving k — i vanish, so equating the above with the known quantity it is equal 
to gives a system of equations (with iu fixed) for the rxy,iu,kw with k ^ i. The jv, kw entry 
of the matrix for the system is 



(-l)"" {Ti - Tk) (jj - Tk) Tkwjv + ijj - Ti) ijv SjV; 



kw 



The determinant of this matrix is a polynomial in the entries of r and ^. (We think of 
the Tj as fixed.) In the expansion of the determinant one summand is Ylj,u{Tj — Ti) ^jy. Every 
other summand will contain at least one Vkwjv factor. If we look at the series expansions for 
these other summands valid for small (coming from the series for the Bessel functions 
and the Neumann series for the resolvent), every term will be a product of powers of the 
and have as coefficient a negative integral power of r(a) times a rational function of a. It 
follows that in the series expansion of the determinant the coefficient of Yljy ^j^, is nonzero. 
Thus the determinant cannot be identically zero. 

We have shown that ifi ^ j then TxySujv is expressible in terms of the unknown functions. 
It remains to consider the cases where i = j, and we always evaluate at {^iu,^iv)- In this 
case (7.16) shows that Rxx — Ryy is known. Subtracting Ti times (7.16) from (7.17) gives 

= xRx + yRy + Rt5Rx - RyTSR - Ti{R5Rx - Ry5R). 

All terms here involving S are sums over k. The terms involving k ^ i, even after taking dx 
or dy, are known, as we have shown. Those involving k — i cancel, just as before. Hence 
applying dx and dy to the above and evaluating at (^j^, ^j^) shows that 

Cm ^xx,iu,iv ~l~ ^iv ^xy,iu,iv and ^xy,iu,iv ~l~ Cif ^yy,iu,iv 

are known. Taking v = u shows that both rxx,iu,iu + rxy,iu,iu and rxy,iu,iu + ryy,iu,iu are known. 
If uj^ V, using the fact that ryy^iu,iv - rxx,iu,iv is known we see also that rxx,iu,iv and rxy,iu,iv 
are individually known. 
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All that we have left to show are (7.6) and (7.10). For these we use (7.15) (the analogue 
here of Lemma 2) and the facts 

(D^ + M-^) = -VP, {D'^ + /5 (1 - /?) M'^) M^j = 2ip - M^j, 

which follow from the differentiation formulas. (The first is just the differential equation 
satisfied by $a; the second is a miracle.) We use these to compute QiujiCiu) and PlujUiu) as 
for previous equations. Thus, for example, to obtain (7.6) we replace the term (^^— 2n— l)g 
in (5.28) by d§ {-p (1-/3) ^'"^q-q) and to obtain (7.10) we replace the term {^'^-2n+l)p 
in (5.32) by d^ (— /3 (1 — /3) C~^P + 2?— p)- Any reader who has come this far can easily supply 
the details. 
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